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We explore the possibility of passive error correction in the toric code model. We first show that 
even coherent dynamics, stemming from spin interactions or the coupling to an external magnetic 
field, lead to logical errors. We then argue that Anderson localization of the defects, arising from 
unavoidable fluctuations of the coupling constants, provides a remedy. This protection is demon- 
strated using general analytical arguments that are complemented with numerical results which 
demonstrate that self-correcting memory can in principle be achieved in the limit of a nonzero 
density of identical defects. 
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Classical information can be reliably stored by encod- 
ing it in long-living metastable states of a many-particle 
system (e.g., the magnetic surface of a disk). The code 
states are typically states carrying different values of an 
order parameter, separated by an energy barrier that 
grows (in the optimal case linearly) with the system size 
[H, Consequently, the probability that an error 
occurs decreases exponentially in N , rendering a high ro- 
bustness of (even very small) memory devices. Crucially, 
the stability of such devices is due to the intrinsic (local) 
interactions between the particles, and no active error 
correction is required during storage. They are thus of- 
ten called 'passive' or 'self-correcting' memories. 

The situation appears far more challenging when it 
comes to the storage of quantum information. While 
quantum information can in principle 47| be protected 
against disturbances using active error-correction 0-0], 
to date no realistic many-body system is known to pas- 
sively (i.e., by virtue of its natural dynamics) preserve 
quantum information over macroscopic timescales. Nev- 
ertheless, self-correcting properties may be obtained by 
a clever design of the Hamiltonian of the system. In a 
pioneering paper @,|3|j Kitaev proposed the toric code as 
a topologically protected quantum memory, where infor- 
mation is stored in degenerate (and locally indistinguish- 
able) ground states. His proposal prompted an intensive 
study of the use and limitations of such systems for in- 
formation storage (see, e.g., Q for an overview). 

In the toric code excitations out of the ground state 
space can be removed by active error correction. How- 
ever, if the corresponding defects have traced out un- 
contractiblc loops when annihilated, this causes an er- 
ror. In particular, it has been demonstrated that Ki- 
taev's 2-dimensional toric code cannot provide robustness 
against the destructive influence of a thermal environ- 
ment or against Hamiltonian perturbations which lead 
to a random walk of the thermally excited defects [l|, 
[lij . Furthermore, these impossibility results have been 
generalized to a wide class of 2-dimensional lattice sys- 
tems 0, On the positive side, a number of vari- 
ants of topologically protected systems have been pro- 



posed where the relevant energy barrier grows with the 
system size 



In this Letter, we analyze passive error correction in 
the presence of a large class of spin interactions or the 
coupling to an external magnetic field, as well as unavoid- 
able or engineered fluctuations in the toric code coupling 
constants. We demonstrate that coherent defect propa- 
gation [i^ would have lead to logical errors, had it not 
been for the Anderson localization induced by the fluctu- 
ations in the coupling constants. In proposed realizations 
of the toric code as an effective model, e.g., from an un- 
derlying Kitaev's honeycomb model |22l |. perturbations 
at the 'physical level' (e.g. dipolar interactions in the 
honeycomb model) induce complicated perturbations at 
the effective level of the toric code, which still can be ap- 
proximated within the class of interactions considered in 
this Letter. The conclusion that the stability of the toric 
code can be improved by randomness was also reached 
in [23| where it was shown that the stability of the topo- 



logical entanglement entropy is enhanced by the presence 
of random magnetic fields. 

Error Model — Consider a square lattice F, with spins 
sitting on its edges, embedded in an arbitrary two- 
dimensional manifold. Assume the dynamics is described 
by the toric code Hamiltonian 



TC 



(1) 



jes 



[Jm > 0, Je > 0, see 01) where the sums run over all 
plaquettes p and over all stars s, respectively. Eigen- 
vectors IV') of Htc with Y\.j(zpZj\'il)) = — IV-") a'l'G said to 
have a 'magnetic defect' at plaquette p. Analogously, an 
'electric defect' at star s is detected by Oj-gs^'^i- 
fects of same type are bosons among themselves. De- 
fects of opposite type have mutual exotic braiding phases. 
This setup serves as a model for topological quantum 
memories. However, it will never be possible to realize 
the iJrc'-dynamics perfectly in any experiment. In par- 
ticular, time-independent perturbations will always be 
present in any experimental realization of the toric code. 
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Two physically meaningful examples of such perturba- 
tions arc the dipole-dipolc interaction 



= 2^ TT^ - 3 rrrji 



(2) 



(rj collects constant factors) with := r.i — fj, and the 
effect of an external homogeneous magnetic field 



H, 



magnetic 
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(3) 



where <7j = (Xj , Yj , Zj ) denotes the vector of Pauli ma- 
trices at spin/edge j and fj collects constant factors. In 
this Letter we consider the effect of translation invariant 
perturbations of the form 



Hr 



E 



C(ii,ai 



,arn^{x ,y ,z} 



(4) 

Here we sum over all collections of m distinct spins / edges 
ii,...,im and all possible choices of corresponding Pauli 
operators ai,...,am G {x,y,z}. The function ^ has 
a spatial cutoff, such that ^(ii, i„i) = whenever 
max{||zfc — : k,l g {!,..., m}} > R {R < oo). The per- 
turbations of the form (jlj can be seen as simple models 
for interactions of the dipolar type, but also include the 
perturbation -ffmagnotic in Eq. ^ . The perturbed Hamil- 
tonian reads H = Htc + Hj. It has been shown [24 1 
that perturbations of the form (|4]) lead to a splitting of 
the ground state degeneracy that is exponentially small 
in the system size. 

In lowest order degenerate perturbation theory one 
approximates the Hamiltonian H by the operator 

J2n,,nr„ Pri^nr^HPn^n^, whcrC Pn,,n^ iS the projcCtOr 

onto the unperturbed eigenspace carrying pairs of 
electric defects and rim pairs of magnetic defects. The 
pairwise orthogonality of the projectors Pn^um allows for 
separate treatments of the dynamics that is induced by 
each of the operators Pn^n^HPn^n^- Consider the sim- 
ple special case rig = 1, rim = 0. To figure out the 
nature of the dynamics described by Pi^qHjPi q we com- 
pute the matrix elements of PifiHjPi Q with respect to 
the basis built up by the toric code eigenstates living in 
the image of Pi.o- Terms that either contain an Xj or a 
Yj = iXjZj give no contribution because they create new 
magnetic defects (or annihilate existing ones) and conse- 
quently map states in the image of Pi^ to its orthogonal 
complement. For the computation of the matrix elements 
we are left with the following expression: 

n n " E £,{ii,z,...,im,z) 

s€li tGb il,---,im'-iq=^it 

X {gi\Zs^ ■ ■ ■ Zs^ Zi^--- Zi^ Zt^--- Zt^ \g.j), (5) 

where 

C := ^ £,{ii,z,...,im,z)Zi^---Zi^. (6) 

ii,...,im:iq^it 



and li = {si,...,Su}, I2 = {ti,---,ty} are paths on the 
lattice, and where the vectors {\gj)}j form a basis of the 
ground state space (the vectors of the form Jlse/a ^slSj) 
build up a basis in the image of Pifi] cf. Q). The matrix 
element within the sum is nonzero if and only if the paths 
li and I2 differ by the movement of the electric defects 
that is determined by the action of Zi^ ■ ■ ■ Zi^^ . Note that 
PifiHiPifi does not allow the particles to sit on top of 
each other (this corresponds to the annihilation of de- 
fects). It thus follows that Pi qHiPi q describes a finite 
range hopping term with an effective hardcore repulsion. 
The evolution of the pair of electric defects determined 
by the spin-Hamiltonian Pi.oHjPi o can equivalently be 
described by a Hamiltonian T for the evolution of the 
defects on the lattice F. These observations generalize — 
up to exotic braiding phases between defects of different 
type — to arbitrarily many defects [4§| . In case of a homo- 
geneous external magnetic field along the z-direction [H^ 
we have an intuitive understanding for the nature of the 
induced dynamics (nearest neighbor hopping of the elec- 
tric defects with hardcore repulsion): the wave function 
of a nearest neighbor pair of electric defects will spread 
arbitrarily during time evolution (cf. the inset of Fig. [1} . 
Note that this leads to the failure of active error correc- 
tion (fusion of nearest neighbors) at the read out, and the 
memory will thus become unstable. Next, we will show 
that this problem will be present for all perturbations of 
the form (|4]) . To apply methods from spectral theory we 
assume F to be Z^. 

Propagation to Infinity — The goal of this section is to 
show that there exist initial 2-defect wave functions (e.g., 
(ne, rim) = (1,0)) with the property that the two defects 
travel arbitrarily far away from each other. This eventu- 
ally leads to an error in the logical qubits. For this we 
need to determine the Hamiltonians that govern their rel- 
ative motion. We follow the approach described in 
and [l^l and find that the relative dynamics is gener- 
ated by a family of hopping Hamiltonians T(fc)o -I- T{k)j 
acting on the relative Hilbert space that is param- 

eterized by the quasi-momentum k {k £ [0,27r)^). The 
operator T(k)Q describes translationally invariant, finite- 
ranged hopping on the relative configuration space 
(i.e., the set of vectors that connect the two defects). 
The interaction T{k)i on the other hand dictates inho- 
mogeneous hopping and is only supported within a finite 
neighborhood of the origin in the configuration space of 
relative motion. In psj we prove that there exist initial 
states such that the two particles are not found within 
finite relative distance as time approaches infinity: 



Theorem 1. For any quasi-momentum k there exist 
e for the relative dynamics 



initial states 



such that lim* 



-i(T[k)a+T[k)i)tj,{k) 







for any finite subset A of the relative configuration space 
1?. 
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Effective Random Potential — Define 



Hr:=-Yl Jrr^ (p) H " E '^^ U ' 



where Jm{p) and Je(s) are positive iid random vari- 
ables described by a bounded and compactly sup- 
ported probability density. In lowest order degenerate 
perturbation theory the total Hamiltonian now reads 

E„„n„ P^.ri^ (^nen™ + ^^r + Hj) P„^„„ . The tCrmS 

HiPn^n^ have been discussed before. We are left 
with the computation of the matrix elements of the op- 



erators Pn 



H P 



with respect to the cigcnbasis of 



the positions of the magnetic and electric defects. Since 
the vectors in this basis are automatically eigenvectors of 

Pn^n^HrPn^n^ WC COncludc that Pn^n^HrPn^n„^ ^CtS aS 

a multiplication operator on the elements of the position- 
eigenbasis. The operators Pn^n^HrPn,,n^ thus play the 
role of potential terms. The potential felt by electric de- 
fects moving on F is specified by iid random variables 
Vs = 2Je(s) (and analogously for the magnetic defects). 

Localization — From the work emanating from Ander- 
son's discovery in 1958 we expect that the iid potential 
that is caused by the randomization of the coupling con- 
stants leads to Anderson localization of the two defects. 
In [2^ Aizenman and Warzel proved dynamical localiza- 
tion of interacting n-body systems {n < 00) on Z,"^ under 
the assumption that the interactions are described by in- 
teraction potentials with finite range: 



(8) 



(x ~ (xi, X2, a;„) G "L ). Here, Aj, V{x), and 
W(x; a) denote the discrete Laplacian, iid ran- 
dom one-particle potential, and interaction po- 
tential, respectively (cf. [29[). They prove dy- 
namical localization with respect to the Hausdorff 
pseudo distance measure defined by dist-H(x, y) := 
max{maxi<i<fc dist(xi, {y}), maxi<i<fe dist({x}, j/i)}. 
As discussed in [2^ this result has the drawback that 
it still allows a particle to hop from one tight cloud of 
particles to another. However, such processes appear 
rather unphysical, and results in [s^, HH suggest that 
they are unlikely. Note that the Hamiltonian ([5]) is 
not exactly of the form ([5]) because in our setting 
the finite-range interactions between the two defects 
are given in terms of inhomogeneous hopping matrix 
elements. To get localization bounds for our 2-particle 
system the proof of Aizenman and Warzel needs to be 
adapted. The details are given in We arrive at the 
following theorem about dynamical localization. 

Theorem 2. Let iJ^^-* he the random Hamiltonian de- 
scribing the evolution of the pair of electric defects on 
the lattice 1? that corresponds to the spin-Hamiltonian 
PiA^i + >^Hr)Pifi. For each m£N (cf (g]) ) there is a 
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FIG. 1: (Color online) Left: Anderson localization of the two 
particles (A — 50, L = 10, tmax = 60). Inset: Dynamics 
of the two particles in absence of a random background po- 
tential. Right: Comparison between the 1-particle dynam- 
ics with and without the exotic braiding phases (A = 50, 
L = 40, tmax = 100). The figure suggests that localization is 
not greatly aflected by exotic braiding phases. 



Ao G ]R.-)_ with the property that for all A > Aq there exist 
A, ^ < 00 such that for all x, y G Z** 



E 



sup |(y|e" 
teR 



< Ae 



-distH(x,y)/5 



(9) 



Neglecting the exotic braiding phases acquired by 
braiding defects of opposite type the theorem general- 
izes to finitely many defects. Taking into account such 
phases is still an open problem. However, if only one 
defect type is dynamic (e.g., in case of @ with B along 
the z-direction) the static defects influence the evolution 
of the dynamic defects in terms of random vector po- 
tentials. First steps towards the proof of localization for 
such systems have been taken in (32| and [33j . 

Numerics: localization of two electric defects — 
Consider a. L x L lattice with periodic boundary con- 
ditions. The dynamics is described by nearest neigh- 
bor hopping together with a background potential that 
is uniformly distributed on [0, A]. The left part of Fig. [1] 
displays the 2-particle dynamics for a typical potential 
landscape. At time t = the 2-particle wave function is 
a position-eigenstate |x), x = {xi,X2) with particles "1" 
and "2" being nearest neighbors. For each y = (2/1,2/2) we 
record the value sup^g^^^ ^^^^^j |(y| exp(— ifft)|x)| and 
plot it as a function of the relative distance II2/1 — 2/2II 1 — 1 
(II • 111 denotes the 1-norm; II2/1 — 2/2II 1 = 1 corresponds 
to nearest neighbor conflgurations) . Each point in the 
figure thus corresponds to a specific 2-particle configura- 
tion (2/1,2/2) with 2/1 and 2/2 on the periodic L x L lattice 
under consideration (L = 10). 

Numerics: influence of the exotic braiding phases — To 
investigate the effect of braiding phases on localization, 
we consider a single electric defect in the presence of two 
static magnetic defects. The pair of magnetic defects 
are generated by a path operator Y[j£i -^j applied to a 
vector in the ground state space, with I being a path 
on the lattice. The nearest neighbor hopping terms that 
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n 



FIG. 2: (Color online) Plot of the critical disorder Ac as a 
function of the particle density n. We assume nearest neigh- 
bor hopping (e.g., caused by (|3])) together with an iid back- 
ground potential whose values are uniformly distributed on 
[0, 2A]. The boundary conditions are periodic and the unit is 
the hopping t. 

cross I change sign. In the right half of Fig. [l]we compare 
the localization in presence and absence of such a path 
along a straight line of length L/2. The localization of 
the electric defect appears unaffected by the presence of 
the path, and thus of the exotic braiding phases. 

Numerics: localization at positive densities — To nu- 
merically estimate the stability of the perturbed toric 
code at positive densities (i.e., infinitely many defects 
in the infinite system) we use Quantum Monte Carlo 



worm-type simulations [3J| (here in the implementation 
of Ref. [35|) to determine the phase transition from the 
superfluid phase to the insulating Bose glass phase (cf. 
Fig. HI, similar as was done in Ref. [s^. The property of 
being insulating suggests that the two defects of a pair 
do not travel arbitrarily far away from each other. The 
details are given in (isj . 

Conclusion: Our results suggest that logical errors 
caused by interaction-induced propagation of defects can 
be suppressed by randomness in the toric code Hamilto- 
nian. While the mathematical approach was applicable 
for a broad class of perturbations and finitely many par- 
ticles in infinite systems. Quantum Monte Carlo simula- 
tions allowed us to draw conclusions for infinite systems 
with positive defect densities that are exposed to a homo- 
geneous magnetic field along the z-axis. Single particle 
simulations indicate the firmness of our findings against 
exotic braiding phases. 

Additional Note — Similar results have been obtained 
independently in [sl] by J. Wootton and J. Pachos. 
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SUPPLEMENTARY MATERIAL 
TORIC CODE 

Consider a square lattice F C embedded in an arbitrary two-dimensional manifold with spins sitting on its edges 
and assume the dynamics being described by the stabilizer Hamiltonian 

p j£p s jes 

{ Jrn > 0, Je > 0) whcrc the sums run over all plaquettcs p and over all stars s respectively. This exactly solvable 
model bears the name Toric Code Its groundstate eigenspace 

Eig{HTc,Eg) span{|.gj)}j=i^...^22s+'» (11) 

is 2^^'+''-times degenerate if the spin-1/2 square lattice is embedded in a manifold with genus g and h holes. The 
operators 

Wlr^ := JJXj (12) 

(products of X- or Z-Pauli operators along the paths / on the square lattice F and I* on the dual lattice F*) are 
sometimes called electric and magnetic path operators. They are convenient to describe excited eigenstates of the 
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Toric Code Hamiltonian Htc because each i^TC-eigenspace is spanned by vectors that result from the apphcation of 
some electric and magnetic path operators to vectors in the ground state space: 



Eig(ifTC, E„^,nJ - span {pf/^' • • • wl' 



open 



paths I 



(13) 



There is a lot of redundancy in this describtion of the eigenspaces because two vectors W^^^ \gj) and Wi^'^ \gj) are equal 
if (1) the endpoints of li and I2 are equal and (2), if the paths li and I2 are homotopic. All electric or magnetic path 
operators along closed paths that can be contracted to one point act as the identity on the ground state space. On 
the other hand, electric and magnetic path operators that cannot be contracted to one point may cause a non-trivial 
linear transformation from the ground state space to itself. The energy-observable Htc is everywhere indifferent to 
the electric and magnetic operator strings except at their endpoints. This leads to the jargon that there are electric 
charges sitting at the endpoints of electric path operators and magnetic charges sitting at the endpoints of magnetic 
path operators. The commutation relations Pauli matrices lead to anyonic statistics between charges of different type. 



ERROR CORRECTION 



The quantum information that needs to be stored with the help of the toric code is encoded in the degenerate 
ground state space. If the spin system lives on the torus the ground state space is 2^-dimensional and can thus carry 
2 qubits. Assume that we have encoded our quantum information into the ground state space at time t ~ 0. As time 
passes the toric code interacts with its environment and the state will acquire support outside of the ground state 
space due to the interactions with the environment and the imperfect experimental realization of the toric code. Error 
correction may move the state back to the ground space. A possible algorithm for error correction goes as follows: 

1. Measure Htc- This is equivalent to the measurement of all the plaquette observables —Jm Yijep ^^'^ 

star observables — Jg Y[j£s -^j ■ There are magnetic charges associated to the plaquettes p where the measurement 
of — Jm rijep yields +Jm and there are electric charges associated to the stars s wherever the measurement 
of ~Je rijgs yields +Je. 

2. Pair up all the charges of same type in such a way that the sum of the relative distances between the 2 charges 
of each pair is minimal. 

3. Fuse the 2 charges associated to each pair by the application of electric path operators Ty/'^'* and magnetic path 
operators Wjl™^ along shortest paths connecting these 2 charges. 

Non-contractable electric and magnetic path operators can form non-trivial maps from the ground state space to 
itself. Therefore, a possible scenario in case of which the above error correction fails goes as follows: consider an 
embedding of the toric code into the surface of a torus and assume that sometimes after the initialization of the 
memory the interactions with the environment causes a pair of magnetic charges. As time passes the imperfect 
experimental realization of the toric code leads to a movement of the particles in such a way that after a certain 
time the particles may have moved more than halfway around the torus. Consequently — after the fusion of the 2 
particles in step 3 of error correction — the ground state vector at t = and the ground state vector after the error 
correction differ by the application of a magnetic path operator along a path that cannot be contracted to a point. 
Such a path operator acts non-trivially on the ground state space and we thus end up with an error in the stored 
quantum information. 



PROPAGATION TO INFINITY 



The goal of this section is to show that there exist initial 2-defect wave function with the property that these two 
defects travel arbitrarily far away from each other during time evolution. This will be achieved by showing that the 
absolutely continuous spectrum of the Hamiltonians generating the relative motion of the two charges is nonempty 
which — via the RAGE theorem — implies the assertion. Let us briefly recall these notions. According to the spectral 
theorem there exists a spectral measure {Pa}acr for every Hamiltonian T such that 



T= f XdPx. 



(14) 
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Moreover there exists a decomposition 

Pa = P^r^ + Pt^ + Pt^ (15) 

with the property that the measure 

/.^A) Pa (16) 

on M is pure-point if € Im(PA^^^''), absolutely continuous if e Im(Pl^°^'') and singular continuous if S Im(P^'"^''). 

Define the subspaces Hpp = lm.(P^^'^), Hac = Im(-PR°^'') and Hsc = Iin(PR^'^^). The restrictions Tpp := TI-h^^, 
Tac ■= T\-fi^^ and Tsc '■= T\-h^^ lead to the definitions of the pure-point part <Tpp{T) ~ (^{Tpp) of the spectrum, the 
absolutely continuous part aac{T) = u{Tac) of the spectrum and singular continuous part <Tsc{T) = a{Tsc) of the 
spectrum with the property 

cT{T)=<jpp{T)yjaac{T)UcTsc{T). (17) 
The RAGE theorem implies the following: Let ?/'o G T-Lac and A C Z'^ finite. Then: 




In other words — assuming that -0 is a 1-particle wave function — the probability for finding the particle within 
any finite region A C Z"* vanishes as time approaches infinity. Hence, if we were able to show that T-Lac 7^ for the 
Hamiltonians generating the relative motion between the two charges we would know that there exist initial 2-defect 
wave functions which will not be found within finite relative distance as time approaches infinity. Consequently, error 
correction will fail as time approaches infinity. The goal of the remainder of this section is the proof that indeed 
T~iac 7^ for the Hamiltonians generating the relative motion between the two charges. 



Relative Dynamics 

To show that the absolutely continuous spectrum of the Hamiltonians that generate the relative motion between 
the two defects is non-empty we first need to determine these Hamiltonians. For systems that live in the continuum 
R'^ there is a single Hamiltonian Hj-d that describes the relative motion of the two particles with respect to the center- 
of-mass frame. On the lattice, the situation is a little more delicate. We follow the approach described in Ref. |27| that 
is itself based on Ref.lH Let (is : l^i'L^) l^i^^^), 

{Usip){ni,n2) := ip{ni - s,n2 - s) (19) 

(s, ni, € Z^) be the representation of the group of simultaneous translations on the 2-particle Hilbert space P(Z^). 
The discrete Fourier transform A : ^^(T^ x T^) x I?), 

kni,n2) = j^ f d^ki£k2i^{h,k2)e-''^-'' (20) 

(T = [0, 27r)) and its inverse V : /^(Z^ x Z^) ^ L^(T^ x T^), 

V'(fci,A:2) = ^ V'(^l,n2)e*"•^ (21) 

are unitary. On i^(T^ x T^), the representation of the group of the simultaneous translations of the two particles 
reads 

(f/,V)(fci,fc2) = (v 0^7,0 AV-) (fci,fc2) = e*^-('=^+'==V(fci,fc2). (22) 

We conclude that 

^slL2(Ffc) = e'" (23) 
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if 

Ffc := {(fci,/c- fci) eT^ xT^lfci eT^}=T'^ (24) 

for all k £ T'^. Therefore, L^(Ffc) is the isotypical component of the l-dimensional irreducible representation of 
the group of simultaneous translations of the two particles with character e** *^, s G Z^. The operator V o T o A is 
decomposable with respect to the fibration 

^ 2 /'*Tp2 



L^(T^ X T^) = / ®L\¥k)dk (25) 
because [V o T o A, Us] = for all s E 1? . Therefore, T itself is decomposable with respect to the fibration 

f{Z^xZ^)^[ ®A{L^{¥k))dk. (26) 
We can express this observation in the form 

r = / S)T{k)dk (27) 

Note that 

a{T) = y a(r(fc)) (28) 

Thus, to prove our claim aac{T) 7^ from above it suffices to show that (Jac{T{k)) ^ for some k e T^. To figure 
out the spectra of the operators T{k) we need to compute the operators T(fc) more explicitly by the determination of 
their matrix elements with respect to a family of wave functions that generate the spaces A (L^(Ffe)) C ^{7? x 7?). 

Lemma 3. Let k e be arbitrary. The set Bk := {^^f ^ J ^ 



6^^7.1,7.2) := (^e— '=^,,(„,_„,) (29) 



(ni,ri2 G li"^ ) forms a complete orthonormal basis in A (-L^(Ffc)) C /^(Z^ x 1?). 

Proof. The family {e*-'i«ie*-'29=|ji, j2 e Z^} of functions (91,92 e T) form the standard basis in L'^{T'^). The map 
Uk : L^(Ffc) ^ L2(t2) defined by 

iuk4'){q):=V2^P{q,k~q) (30) 

(g G T^) is unitary, i.e., 

(■01, V'2)L2(Ft,) = (ufcV'l,Wfc'02)L2(T2). (31) 

Here, 



(V'i,7^2)l2(f,) = / d4k7^i(k)xF,(k)V2(k)xF,(k) (32) 

i[0,27r)i 

and XFfc(') denotes the characteristic function with respect to F^ C T^. Thus {u^^e^^'^^'^e^^^^'^lji, j2 S Z^} is a 
complete orthonormal basis for L'^{¥k). These observations allow for the explicit computation of the wanted basis 

Bk := {(Aow-ie''^i(-)e'-''^(-)) (711,712)! J G Z^}. (33) 

The calculation yields the functions in Eq. (P^ . □ 
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We continue with the computation of the matrix elements 

T{k){l,j) = (ep\T(fc)ef ) = (e|'=\ref ) (34) 
e Z2, fc e T2) of T{k) with respect to the basis Bk in A (^^(Ffe)) C ^(Z^ x Z^). Note that 

supp{ef (•, ■)} = {{71^,712) e X Z\] = 711 - "2)}. (35) 

From the perturbation theory stems the constraint that the particles are not allowed to sit on top of each other 
because this would lead to an annihilation of the two particles resulting in a departure of the state vector from the 
unperturbed eigenspace. We can realize this boundary condition dynamically by setting all matrix elements of T 
that would lead to double occupations of lattice sites to zero. From ([55]) we conclude that the forbidden subspace is 
spanned by {e^^^lk G T^}. To dynamically realize the boundary condition we thus set 

T{k)il,j)^0 (36) 

whenever / = or j = 0. We are thus left with the computation of the matrix elements 

r(A.)(Z,j) = (ep\Tef) (37) 

for I 7^ and j ^ 0. Let Tq be the translationally invariant and self-adjoint operator (arbitrary translations in Z**) on 
P(Z'^ X Z^) with the property 

ro(ni,n2,mi,77i2) = T{ni,n2,mi,m2) (38) 

(tt,!, 71,2, mi, m2 G Z"^) whenever \\ni — n2\\2 > i? + 77i + 1 (i? denotes the spatial cutoff of the spin-spin interaction; 
recall that max{||zfc — 7/|| : k,l ^ {1, 771}} > R implies zero interaction). Note that Tq is well-defined because there 
are no interactions (between the two particles) within the constraining domain 

{(711, 712) e Z-^ I ||7ii - 712II2 > i? + TO + 1}. (39) 

Define 

fr.= T - To. (40) 

We conclude that Tj{-.-) is only supported on a finite neighborhood of the origin in Z"^ and that Tq and f/ are of the 
form 



(roef)(n) = E4°)ef(n-p) 
pev 

(fjef) (n) = E (m - 712) ef^ (n - p) (41) 

(S^^^'' is by construction translationally invariant while ^p^''(7ii — 712) is not). The explicit form of the constraining 
domain ([5^ implies 

C^(71l-712)=0 (42) 

whenever ||7ii — 712II2 > i? + 7n + 1. The set T' C Z^ is contained in {a e Z'*|||a||oo < tii + 1} and the coefficients 
^p*^ + ^f/"* are equal to the sum of those coupling constants ^(ii, z, i™, z) with the property that the difference 
between the endpoint and the start-point of the path that corresponds to {ii, ...,i,„} equals p G Z^. Thus, 

T{k){l,j)={4''\Te^''^ 



I .--J 
ep',(fo+f,)ef ) 



= 5] #^(("1,^.2)) (e (4°^+^p("i-"2)) ef ((711,712) -p) 

\pe-P 



(ni,ri2)GZi 

\pev 



(27r)4 
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The substitutions Ur := ni — n2 and Ha '■= TI2 lead to 



r+P2-Pl 



(2^) 



% I E e'- '^e— ^ I I E (4"' +4(0) e^^'-'=<5,,+,,_p, | (44) 



The divergent first factor is the scattering part that we ignore by setting it equal to 1 because we are interested in the 
relative dynamics of the two particles. This divergent part also appears in the two-body problems in the continuum; 
it corresponds to the free evolution of the center of mass. We continue by reinterpreting the space that carries the 
labels I and j as the configuration space associated to the relative dynamics of the two particles. We do this for each 
fc e separately. The relative dynamics is thus generated by the hopping Hamiltonian r(fc)o + T{k)j where 

T{k)o{l,j) := 7^ E 4°^ e^^^-'-J.v+p.-p, (45) 

T{k)i{l,j) := 7^ E 4'^(0 e^''-''hi+P.-v^ (46) 
^ ""-^ pep 

Note that T{k)j{l,j) vanishes whenever l>R + ni + loij>R + m+l. This is a consequence of (|^^ and the self- 
adjointness of T{k). At infinity the relative dynamics is thus generated by translationally invariant and finite-ranged 
hopping (see 



Existence of Nonempty Absolutely Continuous Spectrum 

In the present subsection we show that the absolutely continuous spectrum aac{T{k)a + T{k)i) of the family of 
operators describing the relative motion of the pair of defects is nonempty. To that purpose we first show that 
(T(T(fc)o) = aac{T{k)Q) and use the Kato Rosenblum Theorem (that is stated below) afterwards to infer 

^ nac{T{k)o) C aac{T{k)o + T{k)i). (47) 

Lemma 4. The spectrum ofT(k)Q is absolutely continuous, i.e., 

a{T{k)o)^aac{T{k)o). (48) 

Proof. Fix k G [0, 27r)^ and set H := T{k)o. According to the spectral theorem there exists a projection valued 
measure (Pa) A such that 



/ rdPr. (49) 
To prove the assertion we have to show that for every i/j e there exists a function e L^(R) with the property 

M^(A) := {^,Pa^) = {i^, jjPrm^ j^m^{r)dr (50) 
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(A C M). To determine dPr{il}) we have a closer look at Hip: 

2 

'[0,27r) 



m){n) = (^i-^ 1^ ^ V(<Z) (i/e-^-? * )) in) dq 

^Piq) (i?,e-^^-(-)) (71) 



1X2 



[0,27r)2 

V'(<?(r,0))ii;,(,.^)e-'^('''^)-"d0rdr 

drr 





^X2 

27r 



{(r»|g(r,0)e[O,27r)2} 



((r, (/)) e R+ X [0,27r) are polar coordinates) with 



^"-(AiT. 4°^ e'^-'=e'^-(^=-^^) (51) 
^^^^ pep 



because 



(r(fc)oe'?-(-))(n) = ^ T(fc)o(n,m) 

2 



pGP 



(27r)4 



gip2-'i:gig-("+P2-pi) 



(27r)4 
^ ' per 



From Eq. (|49p we thus conclude 
(A C R). Now we can go back to Eq. ((50)) : 



dPrW = J^Xio^2V2.]ir) J ^ ^ (i) '/'(9(r,0))i?,(.»e-^«('^'«-"#dr (53) 



(^,y dF,(^)) = ^ V;(n) fj dP,.{^P)\ (n) (54) 
X[o,2V2.]ir) J (^^) V'('?(r,0))£;,(.,0)fE ^(n)e'^«(-.^)-«jd0dr. (55) 



Consequently, 



ii', J^dPrW) ^ J^X[oaV2.]ir) J (^^) i?,(.,^)|V(g(r,0))pd^dr (56) 



'{0(r)} 

because the bracket equals the inverse discrete Fourier transform. The definition 

l{<t>(r)} 

concludes the proof of the Lemma. 



■■= X[o,2V2.]ir) J (^^^ Eg^r,,p)mq{r,^))\'d^eL\R) (57) 



□ 
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One says that the generalized wave operators Q^{A,B) exist if the strong Umits 



n^{A,B)^s- \im e*^*e-*^*Pac(S) (58) 

>-=poo 



exist (see Ref.|38|). The operator A \ Ran 17=*= (A B) is unitarily equivalent to B \ [Pac(-B)'^(Z^)] if ^''^{A, B) exist (see 
the proof of part (c) of Proposition 1 in Ref. |38j) . In our case, the Kato Rosenblum Theorem guarantees the existence 
of the operators r2*(T(fc)o + T{k)j,T{k)Q) (recall that T{k)j is trace class): 



Theorem 5 (Kato Rosenblum, see Ref. |38I ). Let A and B be two self-adjoint operators with A — B being trace class. 
Then, il^ {A, B) exist and are complete. 

We conclude that T(fc)o \ [Pac(r(fc)o)?^(Z2)] = T{k)o (see Lemma g]) is unitarily equivalent to T(fc)o +T{k)i \ 
Ranl7±(T(fc)o + T{k)i,T{k)o) and therefore 

a{T{k)o+T{k)i rRanf)±(T(fc)o+T(fc)/,r(fc)o)) = a.^ (T(fc)o + r(fc), f Ranr!±(T(fc)o + r(fc)/,r(fc)o)) 

= <Jac (r(fc)o) ^ 

(use Lemma 2]). This proves the following Theorem because 

cTac (T(fc)o + T(fc)j \ Ranf7±(r(fc)o + T(fc)/, T(fc)o)) C ,T„,(T(fc)o + T{k)i) . (59) 

Theorem 6. The absolutely continuous spectrum ofTik)^ + T{k)i is nonempty. 



ONE-PARTICLE LOCALIZATION 



From the work emanating from Anderson's discovery from 1958 we know the following: Let Hq be a one-particle 
Hamiltonian on the lattice Z'' with finite-ranged hopping and assume that ■0 is some scattering state of Hq. Assume 
Hq is perturbed by a random potential XV. Then — when turning on A — all the eigenvectors and generalized 
eigenvectors of the unperturbed Hamiltonian that are associated to some energy interval / C M turn (under some 
circumstances) into exponentially decaying eigenvectors i.e., there exist An such that 

|0„(x)P < A.e-^"l^l (60) 

for all eigenvectors (/)„ of Hq + XV with eigenvalue E Cz L This effect is called spectral localization m / C M. Another 
established notion is dynamical localization m / C M: Every initially localized wave function with spectral components 
in / n cr(Ho + XV) will remain exponentially localized at all times. A sufficient condition for dynamical localization is 

< Ae-l^-^l/« (61) 

for some A,£^ G (0,oo) which depend themselves on / {Pi{H) is the spectral projector of H with respect to /). The 
notation E [...] denotes averaging with respect to different realizations of the random potential XV . Note that this is 
exactly what we are looking for: The expected probability for measuring the initial condition (i.e., 5x) at y decays 
— independent of the time t — exponentially in the distance between x and y. From a naive point of view one 
could expect that the notions "dynamical localization" and "spectral localization" are equivalent. In 1995, del Rio et 
al. showed that this is not true in general: while dynamical localization implies spectral localization the opposite is 
sometimes violated [sjj . However in 1994 Aizenman proved a more general version of the following sufficient criterion 
for dynamical localization (cf. Ref. |40|) . 

Theorem 7 (sec Ref. Ilof ). Let H = T + V{x) be a Hamiltonian acting on the Hilbert space ;2(Z'^) with T being a 
finite-range hopping operator and let H^-i, Q C Z"^, be obtained from H by setting to zero all the hopping terms starting 
of and ending at elements outside ofil. Assume that the values V{x), x € Z'', building up the potential are iid random 
variables with a probability measure given in terms of a compactly supported and bounded density p{v). Then, 

¥.[\{5y,{Hn-E)-Hx)Y] <Cse-^\--y\ (62) 

for all E £ I with some ^-independent constants s e (0, 1), /i > 0, C < oo forms a sufficient condition for dynamical 
localization (|6ip . 



sup I ((5- 



Pi{H)5..) 
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The method to prove dynamical localization via the criterium (|62p is sometimes called Fractional Moment Method 
(FMM) . The criterium (j62[) is satisfied for a large class of physical systems in any finite dimension including systems 
with so called "high disorder" : 

Theorem 8 (see Ref. fill ). Let H,Hfi be as above. Then there exists a Xq < oo such that for all X > Xq and all 

energies E 

E [\{5y, {Hn - E)-'S,)['] < ^e-^'l^-^l (63) 
with n > 0, Cs < oo and \...\ denotes (for example) the 1-norm on TL'^ . 

This theorem is a special case of Lemma 3.2 in Ref. |4l| by Aizenman and Molchanov. Other proofs of localization are 
based on the so called multiscale analysis (MSA) approach invented by Frohlich and Spencer [4^ . We conculde that 
if the perturbation W leads to the satisfaction of the criterium (|62|) then the expectation value of the amplitude for 
measuring a particle (that has been initially compactly supported) outside some ball with radius R goes exponentially 
to zero with increasing R. 



MANY-PARTICLE LOCALIZATION 



Note that the 2-particle system with state space Eig(_ffTCj ^-i.c) ^^nd dynamics Pi,oHjPi o is an interacting two- 
particle system so that we are actually not allowed to blindly apply the 1-particle theorems from before. Luckily there 
has been a lot of progress in the investigation of localization in interacting many-body systems. In 2009 Aizenman and 
Warzel (29l | proved dynamical localization of interacting n-body systems (n < oo) on Z"^ under the assumption that 
the interactions are described by interaction potentials with finite range. More precisely, they considered n-particle 
Hamiltonians with up to m-point interactions (m < oo) of the form 

n 

H^"^=Y,[-^j+^V{x,)]+U{K;a) (64) 

(x = {xi,X2, ...,x„) € Z''") where Aj, V{x) and W(x;a) denote discrete Laplacian, iid random one-particle potential 
and interaction potential respectively (cf. Ref. [i^ for the precise specifications) and prove a more general version of 
the following theorem. 

Theorem 9 (see Ref. [29l ). For each n € N and m G {1, ...jjij there is an open set C M_(. x M™ which includes 



which at some ^, ^ < oo and all (A, a) G L^' , and all x, y G Z 



regimes of strong disorder and weak interactions (see Ref. for the precise characterization of these regimes), for 

< ^g-d«t„(x,y)/«^ (65) 



E 



sup|(<5x,e-^*^'"'jy)| 



.tGl 

where the so called Hausdorff pseudo distance is defined by 



dist-u{x,y) := max < max dist{xi, {y}), max dist{{x.} , yi) > . (66) 

I l<i<fc l<i<k J 

As in the 1-particle case (see Theorems [7| and [8|) , Aizenman and Warzel prove Theorem [9] by deriving first the 
many-body analog of the sufficient 1-particle criterium [7] and by showing the validity of the validity of this criterium 
afterwards. Thus, the first step in the proof of Theorem [5] is the generalization of the sufficient criterium described in 
Theorem [7] to interacting many-body systems: 



Theorem 10 (Fractional Moment Criterium, see Ref. 1291). Let i/^"^ be the Hamiltonian (j64p acting on the many-body 



Hilbert space / C and let (T2 C Z'^J be the finite-volume operator that is obtained from H^"''^ by keeping 

all matrix elements that map C'-"-'(ri) to itself unchanged and setting all other matrix elements to zero. Assume that 
the values V{x), x G Z'^, building up the background potential are iid random variables with probability measure given 
in terms of a compactly supported and bounded density p(v). Then the following is a sufficient criterium for dynamical 
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localization (with respect to the Hausdorff pseudo-distance) for energies within the interval I: There exist A,^<oo 
and s €z (0, 1) such that 



sup sup / E 

7CR rtizTrd \1 I 
|/|>1 



dist-ij (x.y) 

dE < Ae « 



(67) 



The existence of exponential bounds for fractional moments of finite-volume Green's functions (cf. (|67p ) is equivalent 
to the existence of exponential bounds for finite-volume eigenfunction corrolators, i.e., 



sup sup E 

/ci 



Q[r^(x,y;/) 



< Ae « . 



Here, 



in 



(68) 



(69) 



In the remainder we will refer to the sufficient criterium from Theorem [TU] in terms of "fraction moment criterium" . 

PROOF OF THE ri-PARTICLE FORMULATION OF THEOREM 2 

Our goal is to use dynamical localization to stop the spreading of defect-wave functions that we have encountered 
at the example of relative 2-defect propagation in the section labelled "Propagation to Infinity". Theorem IH] does not 
immediately lead to dynamical localization because in our setting the finite-range interactions between the two defects 
that arc evolving according to PifiHiPi^ arc given in terms of inhomogeneous hopping matrix elements. Thus, in 
order to get a localization bound for the many-defect system under consideration we need to go through the proof of 
Theorem [S] and adapt it to the setting described in Theorem 2. The following proof of Theorem 2 not only covers the 
2-particle cases that are described in the main text but also general ri-particle cases (n € N). We are allowed to set 
a = (see Eg . ([64])) because we are not dealing with interactions that are specified in terms of interaction potentials. 

In Ref. |29| Theorem 2.1, Lemma 3.1, Theorem 4.1, Theorem 4.2, Lemma 4.3, Lemma 4.4, Theorem 4.5 and 
Lemma 5.1 lead to the proof of the many-body version of the fractional moment criterium (cf. Theorem llOp . To prove 
that (|67p and (j68p still serve as a sufficient criterium for dynamical localization in our setup (cf. Theorem 2) we only 
need to adapt the proof of Lemma 4.6 in Ref. [29l : 

Lemma 11 (Lemma 4.6). Let il C Z'^ andE>0. Then for every x G C^"^^)- 



V(0) niin{l,(n|A|)"i}(| maxsuppHol 'imax--E) 



(70) 



with /imax '■= iiiaxij_..._i^gB_i_j^ij IC(*i7 ^7 *m, z)\ and maxsuppi/o := IJ^. [suppi/o(k, •) — k] denotes an upper bound 
on the ^.-dependent support o/77o(x, •). 

Note that |maxsuppifo| < oo because |maxsupp7?o| is determined by the (zi, i,„)-movements in Z"'' and thus 
I max suppiJo I is bounded by the number of points in the ball Bm (0) with 1-norm radius m located at the origin. 

Proof. In Ref. the authors used the Chebyshev-type inequality 1r\(__b,_b)(x) < e~*^(e*^ -I- e~*^) to upper bound 
the LHS of (|70p by semigroups. To find convenient upper bounds we continue along the lines in section 3 in Ref. |4^ 
using Duhamel's formula 



e-t(.Ho+v) ^ ^-tv ^ / dTe--''i-Ho)e-^'-^'>'^"°+''^ 
Jo 



(71) 



where Hq and V denote deterministic finite ranged hopping and random potential respectively. The iteration of 
Duhamel's formula gives 



-t{Ho+V} 



E 

m>0 



dri • • • dr™ e-^i^(-Fo)e-("^-"^)^(-i?o) • • • {-Ho)e- 



it-rm)V 



(72) 



0<Ti<...<r„<t 
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As in Rcf. 43 we can rewrite (l72l) so that 



(73) 



Here, V{n) denotes the potential on the n-particle configuration space Z""* that is induced from the 1-particle potential 
XV {u), u £ Ii'^. Note that only the exponential contains random potentials so that 



E 



(^x,e-*(^°+^'<5x)l < / diy{n{T)) 



cxp 



V{n{T))dT 



n 

i=l 



\Ho{xt,yt)\Jn^,{Ti+)ny^{Ti~) 



(74) 



(l-ffol denotes the operator that emerges from Hq when replacing all the matrix elements of Hq by their absolute 
value). Observe that 



E 



exp 



V{n{T))dT 



where 



{s))ds 



(75) 



(76) 



can be regarded as an averaging operation. Pulling it outside of the exponential function and the expectation value 
by using Jensen's inequality we arrive at 



E 



-nt\\\V{u) 



Nu{n{s))ds 



EL-nt|A|V(0)l /V„(n(.))d 



E 



-nt\\\V{0) 



(77) 

(78) 
(79) 



This yields 



E 



< E 

=:E 

< E 



-nt\\\V{0) 



^-nt\\\V{a) 



,-nt\\\V{0) 



di^{n{T))Y[ \Ho{xt,y^)\^Jn^i{n+)ny.{Ti-) 

i=i 



— 1| max supp//o |/in 



(80) 



where maxsuppi^o := Uk I^^'^PP^oC^' ') ^ 1^] denotes an upper bound on the k-dcpendcnt support of iJo(k, ■) 
and ft-max maxn^m l-ff (n, m)| and . A similar bound holds for t > 0. Putting everything together with 
t := — min{l, (n|A|)^^} we get the desired upper bound. 

□ 

To prove dynamical localization of the toric code defects we have to show that the sufficient criterium ([67]) holds true 
for our n-body Hamiltonian (see Theorem 2) which specifies interactions between the defects in terms of inhomoge- 
neous, finite-range hopping matrix elements. To that purpose we need to restate and prove Lemma 5.1, Theorem 5.3, 
Theorem 6.1, Lemma 6.3. and the induction step in Sect. 6 of Ref. [i^. We start with the formulation and the proof 
of the main theorem. Lemma 5.1, Theorem 5.3, Theorem 6.1 and Lemma 6.3 follow afterwards. The main theorem 
is formulated for the presence of electric charges and absence of magnetic charges. Its adaption to the presence of 
magnetic charges and absence of electric charges is immediate. Note that the Hamiltonian Pn^,o(,Hi + XHr)Pn^,o from 
the main text (cf. (4) and (7) in the main text) induces quantum dynamics for the 2ne electric charges on the lattice 
Z'' {d = 2 in case of the toric code). We denote the Hamiltonian that describes this evolution of the 2ne interacting 
electric charges on the lattice Z'' by TJ''^"^'. 
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Theorem 12 (Dynamical Localization). Let n G 2N (n < oo) denote the number of electric charges and let i?^"^ 
be the Hamiltonian describing the evolution of the n electric charges on the lattice Tf' that corresponds to the spin- 
Hamiltonian Pn/2,o{Hi + XHr)Pn/2,o (^f- (4) ^-''^d (7) in the main text). Then for each m G N (cf. (4) in the main 
text) there is a Xq £ (0,oo) with the property that for all A > Aq there are positive constants < oo such that 



E 



.ten 



< ^g-&st„(x,y)/e^ (81) 



for all X, y G (Z^)^. Here E[...] denotes the expectation value with respect to the random potential variables V{x), 
X e Z'^. 



Proof. The theorem is proven inductively in the particle number n; the validity of the fractional moment criterium (|67l 



in the 1-particle case (n = 1, the induction anchor) has been proven in Ref. |41| . Consequently, it is left to show that 
the validity of the fractional moment criterium (|67p for ?i — 1 particles implies the validity of the fractional moment 
criterium for n particles. Without loss of generality we can assume that n C Z'* is chosen such that x, y e C^") (Q) 
(otherwise, the LHS of (|67p vanishes; recall the definition of Hq from Theorem [T0| . Choose Lq G N arbitrarily 
(we will fix the value of Lq at the end of the proof) and set Lk+i := 2{Lk + 1) as in Ref. [i^ Let x,y £ be 
those vectors that realize the HausdorfF pseudo distance, i.e., dist-H(x, y) = |a- — and assume without loss of 
generality that \x ~ y\ > Lq . Consequently, there exists a unique fc G Nq such that y ^ Kl^{x) but y G Ai,^^j(x). 
Assume that x and y arc sufficiently far apart to guarantee that Lk > 4m. Note that there exists c < oo such that 
Lk < dist-j^(x, y) ~ \x — y\ < cLk. The remainder of the proof is divided into two parts. In the first part we assume 
diam(x) > In the second part of the proof wc will thus have to deal with n-defect configurations x G Z"'' with 
the property diam(x) < In case of case diam(x) > Eq. (A. 2) in Ref. [29l i.e. 



yields the lower bound 



l{x) > diam(x), (82) 



;/ ^ ^ ^ A- f \ ^ dist«(x,y) 

((x > -diam(x > — — > — — . (83 



for Z(x). We conclude that in case of diam(x) > the Theorem is a direct consequence of Theorem [T4l because for 
instance in case of Z(x) > l{y) and dist>^(x, y) > Z(x) we get 



sup sup JH/ 
|/|>i " 



G["^(x,y)|" < Ae-t'W < ^ e" ^ (84) 



This observation implies dynamical localization (see Theorem fTU)) . Wc will thus assume diam(x) < in the remainder 
of this proof. Hence, x G C'^"-'(f2 n Al^{x)). With the resolvent identity wc can remove all the terms in Hq which 
connect C^") (17 n Al, (a;)) to its complement in C(")(17), i.e., C^") (17, 17\Al, (x)): 

|Go(x,y,^)| ^\{6^,{Hn-z)-'6y)\ 

- X] K'^x, (-ffc<">(onAi,j^(x)) © ^c(")(a,n\Ai,^(x)) ~ ^)""^'^w) 
w,w'ea(">)Ai,j. {x) 

X (^w,r^;:;j^;yf;/^(^^2^ (85) 

with 

a('")M := |r = (ri,...,r„) G C^^^M) 3s G Z'^'^^C^") (M) such that max{||sj - r^Wi <m\i,j e {l,...,n}}| 

U {r = (ri,...,r„) G Z("'')\C(") (Af) 3s G C^"\M) such that max{||s, - ||i <m\i,j G {l,...,n}}}. 

(86) 
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Here r 
used 



C'">(f2nAi,j£!;)) 



denotes the "boundary strip operator" - -ffc(")(nnAi, {x}) © -ffc(")(f^,SAAi, (x)) and we have 



(q,q')eZ"''xZ"'^|((5q,r 



7id I / 



;:.)^q') ^ C 9(-)Ai,(x) X S(™)Ai,(:r). 



Together with |((5q,r| 



c'-"\nnALjx)) 



q'-^c(")(a,n\Aijx))'- 



^q')| < ^max for all q, q' e Z" we thus arrive at 



(87) 



w,w'e9('"'Ai,^ (an) 

and therefore (recaU that x e n Al^(.t))) 



|Gn(x,y,z)| <Vax ^ l('5x,(-ffc( 



irga('")Ai. Jx)nC<"'(A i,, (x)) 
w'(Ea("«)Ai,^(x) 

E 



I GnnAi . (x) (x, w) I I Gn (w' , y) | . 



(89) 



wga('")Ai. Jx)nC'">(A i.fe (an)) 
w'Ga('")Ai, (x) 



Hence, Theorem 2.1 (the Wegner-type estimate) of Aizenman and Warzel's paper [23| gives (using la + fej" < |a|* + |6|*, 

se (0,1)) 



c 



Ej [\Gn{^,y)n < K...j^ 9(™)Ai,(x) 



E 



Er 



IGjinAz,. (x)(x,w)|' 



(90) 



^ea('")At. (a;)nC(")(At. (a:)) 



Note that 9("')Alj. (a;) grows only polynomially in Lk. Therefore, this factor does not prevent exponential bounds 
in terms of L^- In analogy to p89p we can rewrite the sum to get (see also (|149p ) 

C 



E/[|Go(x,y)r] <h:^,^— 

C 



s(")Ai,(x) J2 E 

«Ga<'")ALfc(x) wGC(")(Ai,^ {x)na;u) 



|GnnA,,,(x)(x,w)|' 



S(HAi,(x) 



Consequently, 



with 



.s^P E E E,[|G,-,(x,w) 

n:i2CAi,,. (a;) „gg(™) (^) (fj.,^) 



G 



supE,[|Gn(x,y)r]</i^^,_ 5(™)Ai,(x) (5i + ^2) 



(91) 



(92) 



|/|>i 



5i:=sup sup E lE,[|G^(x,w)|^ 

|if>i neac") A^ jx) wec(") (n;u)\c<2'^ (s=2;«) 

52 1- sup sup E E/[|Go(x,w)r]. 



(93) 



|GnnAi,,(a;)(x,w)p 



We continue with the estimation of the term Si: By Theorem 1141 

E E 

«6a(") A^,^ (x) weC<") (0;tt)\C^'2'^ (f2;u) 
<; ^ g- J min{dist-H(x,w),max{i(x),i(w)}} 

ti6a('")Ai,j^(x) wec(")(0;ti)\c<"' (n;«) 



(94) 
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with 



dist^(x, w) ;=max< max dist(a;i, w), max dist(x, m^) 

\l<i<k l<i<k 

>dist(x, u) 

>Af^--) (95) 
because u € d^"^^ Al^{x) and diam(x) < by assumption. On the other hand w e C("' {Cl; {Q; u) impUes 



■ n - 1 V 2 



diam(w) > and therefore 



ji-12 ~n-lV2 

The use of ^ and (HI]) in O yields 



Z(w) > >J—(^±- rn ] . (96) 



E E E,[|GonA.j.)(x,w)|-^] < E E Ae-t-^(--™) 



1 1 



<Ae T (97) 

for appropriate definitions of A and <^ because the quantities |(?^™^Aij.(a;)| and jC'-"-' (Ai^)| grow only polynomially in 
Lfc- The assumptions dist-H(x, y) = \x — y\, y ^ Kl^.{x) but y £ Kl^^^{x) with L^+i := 2{Lk + 1) at the beginning of 
the proof imply |a; — y| < 2(Lfc + 1) and thus dist>^(x, y) < 2(Lfc + 1). We conclude that after the redefinitions of ^ 
and A we arrive at 

~ _ dist-^ (x,y) 

Si <Ae e . (98) 
To finish the proof of the theorem we still have to upper bound the summand S2- we get 
52 = sup sup E E E,[|Gf-,(x,w)n 



< sup sup |9(™)AiJ E E |4r (^,y) 



=i?(")(ifc) (99) 

(see (|148p ) using the translation invariance of the expectation values for the first equality. The second inequality is 
simply a consequence of adding a factor > 1 and new summands to the sum. Therefore (use ([^5)) and ([M)) in ([M])), 

Ei [|Go(x,y)r] < h:^,^^ (^Ae'"-''^ + s(«)(i,)^ . (100) 

It is thus left to show that B^J^\Lk) decays exponentially in the distance dist>i(x, y). Because of the assumed uniform 
{n — 1) particle localization Theorem 1151 implies that there exist a, A^p < oo and v > (these quantities depend on 
the localization properties of the n — 1 particle system) such that 

i3i"'(ife+i) < Y^Bt'Huf + ALf_^,e-^'^^K (101) 
To prove the exponential decay of B^^^\Lk) we need the Lemma 6.2 of Aizenman and Warzel's paper [29| : 
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Lemma 13 (Lemma 6.2). Assume that S : R — >■ M, q,b,p,'r] S [0,cxd) and i^, £o G (0, cxd) satisfy 



5(2^£o) < qS{2''-^Cof + b{2''-^CofPe~^''^^''''^"\ (102) 
if > qb + V^, (103) 
1 > qS{Co) + Jy-Cge-"^" =: e~^^« (104) 



S{2''Co) < ie"^^'"^" (105) 

q 



for all fe e No- Then, 

for all fc G No . 
Proof. Define 

RiL) := qS{L) + r^L-Pe-"^ 

and observe that 

i?(2'^/:o) = qS{2''Cn) + 77(2'=£o)''e-"2''^'' 

< g (g5(2'=-i/:o)' + fe(2'"'/:o)'^e-2-(2'=-^£o)^ + ,,2^ (2'=-i£o)^e-''2'^° 



Consequently, 



= {qS(2'-'Co)f + ^qb + n^^J^^ {2'''C,fP (e-^^-^^° 

< {qS{2^-^Ca)Y + if{2''-^CofP (e-''2'^'-'^")' 

< {qS{2^-'C,)f + (r7(2'^-i£o)^e-''2'^'"'^«)' 

< (g5(2^-i/:o) + /7(2''"'/:o)''e-''2-"i£o^' 

= i?(2'^-i£o)^- (106) 

5(2'=/:o) < -^R{2''Cn) < ^R{2''-^Cof < ... 

< -RiCof 
q 

q 

for i?(/:o) = e-^^o. □ 

We proceed by using Lemma [T3l to find an exponential upper bound for Bi^\Lk) {S{Lk) := Bi"'\Lk), Lk :— 
2''{La + 2) as in Ref. l29l ) with respect to Lk. If we define £o := -^o + 2 we get 

= SiLk) = 5(2^ (Lo + 2)) = 5(2'^Xo). (108) 

Recall that we have not fixed the value of Lq so far. Hence, our next goal is to fix Lq such that Cq = Lo + 2 and S{-) 
satisfy the conditions (fTOS]) . (fTOS]) and (fT04|) in Lemma [T3l We start with the condition (fT02|) : 

5(2'^Xo) = 5(2 • 2'=-! (io + 2)) = 5(2Lfc_i) 
= 5(Zfe)=i?(")(ifc) 

A 

= ^5(^_i)2+iLf_,e-2-^^-^ 



a 



= — 5(2*'--^(Lo + 2)y + A{2^-'{Lo + 2)) 
A** 

= ^5(2'=-i/:o)' + i(2^-i/:o)'^e"'"''"^« (110) 
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(note that Lk = ~ 2 for all k). In inequality (|109p we have used Theorem [151 Consequently, condition (|102p is 
fulfilled independently of the choice of Lq. To satisfy the conditions (|103p and (|104p we choose £oi such that 

2 ^ ^ 

(note that a choice of £o ^-Iso fixes Lq because £o = -^o + 2 by definition) and 77 such that 

^ < < :j7pe'' " (112) 
(the existence of this regime is guaranteed by (jllip ). The left inequality leads to the observation that 

fi<? (113) 

is sufficient for the assumption (|103p in Lemma [13] (using the correspondences ^ O q and A ^ b) because (|113p plus 
^ times the left inequality of pi2p gives 

We note that the condition (|113p can be satisfied by choosing A large enough. The assumption (|104p of Lemma [T51 
reads 

^5(/:o)+'7^ge--'^" <1. (115) 
According to the right inequality in (jll2p we can rewrite the LHS of this condition, 

^SiCo) + ^Ule-^^'^ = (£0 - 2) + r//:j;e-''^° (116) 

<^B(")(/:o-2) + i, (117) 

^B(")(£o-2)<i. (118) 

This demand can again be satisfied by choosing A large enough because the quantity si"-* (£0 ~ 2) can be bounded by 

B("n'Co-2) < -^(£0 -2)2'*("~i) (119) 
A* 

(this is a consequence of the Wegner estimate; cf. (2.2) in Ref. [jgh where C is A-independent. This concludes the 
verification of the assumptions of Lemma 1131 Its application yields 

5(2'= £0) < —e"^''^". (120) 
a 

Therefore, 

{Lu) = S(Lu) = 5(2^- (io + 2)) = 5(2'=£o) 

< ^e-^^'-Co ^ ^e-^^'' = ^£-'^(^'=+2)^ (121) 
a a a 

Recall that 

^i^hli^ < L, (122) 
c 

because of our observation that Lk < dist-H(x,y) ^ \x — y\ < cLk (c < 00) at the beginning of the proof. We thus 
arrive at the desired exponential decay 

5^'(ifc) < (~^"^^) e-'^'^'''*«('''y). (123) 

We have thus shown the validity of the fractional moment criterium (|67| for n defects (cf. (jlOOp ) based on the 
validity of the fractional moment criterium for n — 1 particles (this was the induction assumption) . This concludes 
the inductive proof of the theorem. □ 



and arrive at the sufficient condition 
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Theorem 14 (Theorem 5.3). Assume a system with n — 1 particles (n > 2) satisfies the fractional moment cri- 
terium (j67[) for all A € (Ao,oo). Then there exist some s e (0, 1), A,£^ < oo such that 



sup sup E, |G^"^(x,y)P 

/CM ncz<i 
|/|>i ~ 



< ^ g-|- min{(iistH(x,y),max{i(x),i(y)}} 



(124) 



for all A G (Aq, oo) and aZI x, y G 



ydn 



Proof. As in Ref. |29| we assume without loss of gcncrahty that Z(x) > l{y) and set J,K C {1, n}, JU/l' = {1, n} 
such that 



Ux) = mill |.Tj - .Tfcl 



(125) 



Corresponding to the partition {1, n} = JUX there is a division of the total n-particle system into two subsystems. 



We set h'^'^^ = h'^^ + H^' where H^^' and H^' are chosen such that the dynamics of H^'"^' and Hq, 

agree within the subsystems. Thus Hq'^^ emerges from Hci by removing all the interactions between the subsystems 
associated to J and K. With the use of the triangle inequality \a + b]'^ < |a|'* + {bl" Aizenman and Warzel split up 
the expression we want to bound; 



AJ) 



AJ-K) 



]E,[|Go(x,y)r]<]E, |G[/'^^(x,y 



-Er 



G[/'^)(x,y)-Go(x,y) 



(126) 



With Theorem 5.2 in Ref. |29| we get the desired bound for Er |Gjj' ^(xjy)]** using dist^' ' (x, y) > dist-^ (x, y): 



E. 



|G[/'"^(x,y) 



< ^g-distH(x,y)/4 



(127) 



Thus, we are left with the task to bound E/ [lAj**] for 

A:=G[/'^)(x,y,^)-Gj^(x,y,z). 
The next step in the original paper is the application of the Holder inequality: 



E/flA^l < E 



.3e(l + .) 
A 2(l + 2s) 



3s(l + s) r , -1 \ 

< c|Apw^ (E/ IAI'*/^ \ 
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(128) 

(129) 
(130) 



where Aizenman and Warzel used their Theorem 2.1 in the second inequality together with the definition /3 := 2(2+ s) • 
The application of the resolvent identity yields 



A = E G[/'^''(x,u,z)(<5,, (i/o - i7(/'^V^)Go(w,y,z). 

u,wec<")(n) 



(131) 



We thus arrive at 



li [|A|-^] < clAp^WE/ 



<c|A| -wi+w^] 



^ g[/^^)(x,u,z)(5u, (i/o - i/,^/'''^)5w)Ga(w,y,z) 
u,wec(")(o) 



s/2' 



3£(l+£) 

<c|A| 2(2+.) 



^ G[/'^^)(x, u, z)(5„, (if,, - //</'^'Vv 

U 

^GW^)(x,u,z)(J„,(iJ,,-<^^V.. 



,/2 



|Go(w,y,z) 



K/[|Gn(w,y,z; 



2/3 



2/9 



(132) 

(133) 
(134) 
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The last inequality follows from the Cauchy-Schwarz inequality. Theorem 2.1 from the original paper [29[ gives 



Ei[\Gn{w,y,z)\r < G 



(with constants Cs,K,Eo) so that 



C_ 



KEq 



(the constant C collects constant factors). Next we note that ((5u, (i^n — iJ^^'^^')(5w) is not a random variable. 
Therefore, 



(135) 



(136) 



' ' w 

c 



< 



|A| 
CA 



E 



U 



Er 



G'[/'"^^(x,u,z) 



■^\-fdi<'-)(x,u) 



<^ E E|('5u,(i^o-4'''^^^)<5w) 



|A| 



g-fdist„ '(x,u) 



(137) 

(138) 
(139) 
(140) 



uesupp{_ffs2-H*/'^*} 



where we have used Theorem 5.2 in Ref. [H and \a + bl^^ < \a\^ + jS e (0, 1). The hopping we consider is 
finite-ranged. Set m S N equal to the maximal hopping range and let _Bm(u) denote the ball around u with 1-norm 
radius m. Hence, for a fixed u S C'^"^(r2) there are at most |i3m(u)| nonzero matrix elements of the operator Hfi and 
thus of Hq - H^'^\ Therefore, we observe that the bound |i3„i(u)| < (2to)"'' leads to 



CA. 



E.OAII < j^(/i„,ax)^^(2m)"^ E e-T'^'^'n 

u6supp{fffi-//<,^*'>} 

with ft-max := maxn.m |^^(n, m)|. From Eq. (5.19) in the original paper [2^ we know that 



dist^'-^''(x, u) - (/(x) - to) > 0. 



uGsupp{i/r2-_f/<,'-'^'} 



(141) 



(142) 



It follows that 



E 



-|dist<,''^''(x,u) ^ g-|(i(x)-m) 



E 



'■ (dist<^^"^^' (x,u)-(i(x)-m)) 



uesupp{Ht2-_ff</'^'} 



uGsupp{ffn-if^/''^'} 



< g-|(i(x)-m) 



E 



.{J,K) , , 

' ^ ( dist^ (x.u) — mill 



(J.K), dist^'''^'(x,u) 



uesupp{fffi-//^/'^'} 



< g-f (i(x)-m) ^ 

uec(")(f^) 



g-|-dist!^' '(a,u) 



(143) 



where a S Z'"^ realizes the minimum min „(j/f)i distlf'^'^Vx, u). The use of the relation distlf'^''(x, y) > 

dist-H(x,y) yields 



g-|dist(,^-^)(x,u) < g-f(i(x)-™) g-fdistK(a,u) 



uGsupp{//f,-ff^'''^'} 



uec(")(o) 



/3 



(144) 
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The last inequality is a consequence of Lemma A. 3 in Ref. |29i After this little detour we conclude 

Ej[\An<^e-P^-^ (145) 

if we set 



B := CA{K,,^yP{2mr''C{n, d) ( | ) e?™ < oo 



nd 



(146) 



Putting together the Eq.s ()126p . (|127|) and (|145p we arrive at 

Ej [|Go(x, y)!"] < ^e-t (147) 
This proves the theorem because we assumed without loss of generality l(x) > l{y) at the beginning of the proof. □ 
Set 



Er 



with 



i?(")(L,):=|a(™)Ai,| sup sup E 



:= {w e M| min \\q - w\\i < m - 1} 
qed'--Hi 



(148) 



(149) 



where 9^ 'M denotes the inner boundary of M C Z''. Further we define the sequence {Lk}k>o recursively: L^+i :— 
2{Lk + 1). The starting point Lq will be specified later. 

Theorem 15 (Theorem 6.1). Assume that the {n — l)-defect system satisfies the fractional moment criterium (j67p 
for all A G (AojOo), Ao < oo. Then there exist some s G (0, 1), a,A,p < oo, and ly > such that 

a 

for all A G (AqjOo) and all k Cz No- 
Proof Set 



(150) 



^(")(L,+i) |a('")Ai,^J sup ^ sup ^ J2 



Er 



|^|if/2CAi,^, j;ea(")A^,^, x6C<2'^_^^jn;0) 

We follow the strategy of Aizenman and Warzel and write Bi^^L^+i) in the form 
?(«) 



l4"^(x,y) 



(151) 



(152) 



and estimate Bl (Lk+i) and the difference appearing on the RHS separately. We start with the estimation of 



Bi''\Lk+i) - BriLk+i) 



(n). 



i?(")(Lfc+i)-B(")(Lfc+i) = |a(")Ai,^J sup sup ^ 



(153) 



with 
A ~ 



E 



Er 



l4"^(x,y) 



*eC(-> (f2;0) 



E 



l4"^(x,y) 



fc + 1 



E 



(O;0) 



E 



Er 



E i^H...] 



E 



E 



(154) 
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The quantity A can be bounded as follows: 



yeC<"' (f2;y)\C("' (n;y) 

< Si +S2 



E 



xec<';> (n;0)\c<"' (O;o) 
yec'"' (a;y) 



with 



sup 



E 



^2 := sup 

k-y|>2ifc + i-m 



E 



Er 



(155) 



(156) 



(157) 



y6C("' (n;y) 
:eC<"'(n;x)\C^"' _ (0;a;) 



The last inequality a consequence of y € O'^'^^Al^^-^. We proceed by the treatment of ^i. According to (A. 8) in Ref.[2; 



dist«(x,y) > Ice - y| - {Lk - ^m) > 



71-1 



because 



\x^y\> 2Lk+i -m = 2^{Lk + 1) - m. 



(158) 



(159) 



At the same time, Lk/{n— 1) serves as a lower bound for Z(x), so that dist?^(x,y) and max{Z(x), l(y)} are both lower 
bounded by Lk/{n — 1). Theorem [Til thus yields 



|G|"'(x,y)|-' 



< Ae 



(160) 



for the summands appearing in the definition of ^1. The number of summands can be bounded with the help of the 
observation that 



The combination of all these estimation gives 
Similarly, we get 
Note that for example 



d(n~l) ( Lk+l 



d(n-l) 



<((2L,+i)^)^"-^^ = |A,,,J 



We conclude that 



after convenient definitions for A and ^. Going back to (|153p . we find 



i3("'(ifc+i)-s("'(ifc+i) < \^^^-^^L,,^^ A{Ltl-,')\-T^^ 



(161) 



(162) 



(163) 



(164) 



(165) 



(166) 
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after the convenient redefinition of A. One of the factors jc?^™^ Alj.^J emerges from the replacement of the sum 
X]j,ea(")At. (|153|) . To conclude the proof of the theorem wc need to estimate Bi"\Lk+i) (cf. (|152p ). The 

application of Lemma [16] yields the upper bound 



Bi"Hife+i)<7^/^i\x|maxsuppi/on5('"^Ai,^J sup sup la^Ai.^J 

l-f|>i ~ 



(167) 



X sup 



E E E 



yGa(")Ai,^^j yxGCi.2',(n;0) ue9('")(Ai^nO) wGC(")(Ai,^ni2;«) 



g(;^(x,w)|^ 



EE E 

yyec*'2'^ (f2:y) «ea("') (At^^nO) v'6C(") (Ai,^^ (y)nn;«) 



with 1/ = Al^ n fi and = ^L^iy) ^ fi. Wc continue with 

Bi^'HLk+i) <^hllJma^snppHo\^\d^^^AL,^A^ sup | 

i/r>i 

X sup <^ E E 

^^^^fc •ue9(")(Ai nf2) xeC<"' (O;0) 

weC'"'(Ai, nf2:«) 



(168) 



(169) 



i4rHx,w)i^ 



X sup sup < ^ ^ E/ 

v'eC("'(Ai,Jy)nSl;ii) 



|g'"'(x,w)|-^ 



(170) 
(171) 

(172) 



The translation invariancc (simultaneous translation of all particles) allows to shift y E 'E'^ to the origin. Therefore, 



Bi"\L,+,) <^ I maxsuppi7o|\'|a('")ALj2 



X sup sup \] 2_] 

\ /CB OCA, , ~ - 

|/r>i uea('">(Ai,fcnn) xec<2' i^-.o) 

wec<")(Ai nn;tt) 



4"^(x,w) 



(173) 
(174) 



where 77 < oo is defined such that 



(175) 



for all k. We can rewrite this bound in terms of the function fcs(Aij, , rfc, r^) (its definition can be found in Corollary 
5.4 of the original paper): 



SfHifc+i)<^'^C\JmaxsuppiJop(i3i"Hifc) + fc,(A^,,rfc,r,)|aM(Ai,nf2)p) (176) 

(177) 

where m denotes the hopping range. The application of the Corollary 5.4 thus gives 

(178) 



<v'i^hll^ I maxsuppi/oP (si"'(ife) + fc.(Ai, , r^, r,)(m|9Ai J)^) ' 



Bi-HLk+i) < C^ax I maxsuppifoP (s^Hifc) + A^L^^^-^^- 



After the appropriate redefinitions of the constants the combination of the bounds (|166p and (|178l) yields the upper 
bound of the theorem. □ 
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Lemma 16 (Lemma 6.3). Let ft C Z'^ and V,W € n with dist{V, W) > 2. Then for all x G C^"^^) f^^d y G C(")(W^) 
E[|G[,")(x,y;z)r] < ^ /.^^^^ J max suppi?„ p ^ ^ E [jG^^^Cx, w; z)^" 

X 5] ^ E[|G|;^)(v,y;z)r] (179) 

-ugOCmjW vGC(")(W;-u) 

where the constant C — G(s,d) < oo is independent of (A, a) G (cf. (|149p /or f/ie definition of d^^"'\- ■■))■ 

Proof. To remove all the terms in Hq, that connect sites in G'"'(y) with sites in the complement of G'"-'(T^) (i.e., 
G'"-'(ri; ri\l^)) and the analog for V ^ W wc define the boundary-strip operators 

(180) 
(181) 



^ Q.\V 
^ Q.\W 



By assumption X G G(")(y),y G G(")(V^^) and V^nVP" = 0. Note that for x G C^''\V) andy G C(")(W^) C C^''\n]Vl\V) 
we have for example 



{5^,{Hv®Hn\v~z)-^5y)={). 



(182) 



Therefore, the resolvent identity yields 
Go(x,y;z) ={5^,{Hn - z)-Hy) 



= (<5x, -{Hv © Hn\v - z)-^T\y{Hn - zyH^) 

--- {S^,iHv(BHn\v-z)-'S^){S^,TY,\yS^') 



(183) 



w6C<">(y) vec<"'(i^) 
v'ec<"'(a,f2\v) v'GC<"'(n,a\w) 



X ((5w', (-ffn - z)-'S^){S^,T^^^S^,){5^,, {Hw ® - z)-^Sy) 

= - E E Gl7)(x,w;z)rK\^(w,w')G[,"\w',v;z)rjr\M.(v,v')G|:;^)(v',y;z) (184) 

w'GC<"'(n,f2\v) v'ec<")(n,a\H') 
Using Theorem 2.1 of the original paper we thus get 



E[|Gn(x,y;z)r]< E 



E \Tr,\v{y^,^'WK\w{^y)\' 

w'ec<"'(n,o\y) v'GC('''(f!,f2\w^) 
X E [|g(;') (x, w; z)G[,") (w', v; z)G^^^ (v', y; z) 

<7^ E E |rV(w,w')nrjr^^(v,v')r 

wgc'^H^) vec<"'(iy) 
w'GC('''(n,n\i/) v'GC(''>(n,n\iv) 



X E 



G^'^(x,w;z) E G^^(v',y;z) 



We observe that 



{(q,q') G Z"'* X Z"'*|(,5q,rfj\j,fV) ^ 0} ^ X ^^"^^^ 



(185) 



(186) 



(see ([86])). Therefore, rjf^^^( w, w') 7^ for w G n G'-"^(T^) implies w' G d^™''V\C^"'^ {V) , and rj]^^y(v, v') 7^ 

for v' G 9(™)W^nG(")(VF) implies v G a('»)iy\G(")(M^). We conclude that 

G 



;[|Go(x,y;z)n < 



E 



E |rj;^v-Kw')nr,T\w(v,v') 



vGa(")vnc<"'(y) vGa('")w\c<"'(w) 



X E 



G^')(x,w;z) " E G|;;^(v',y;z) 



in), 
W 



(187) 
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Note that the product of the two expectation values is independent of w' and v. Thus, for each w and v' wc have at 
most I max suppffo I no n- vanishing terms that are upper bounded by /imax so that 



C 

E [|Go(x, y; z)|"] <—h'^^^J maxsuppiJo 



E E 



E 



GV"^(x,w;z) E G^^(v',y;z) 



(188) 



Next, we rewrite the sums in the form 



E 



v'ea<'")wnc<")(w) 
to conclude the proof of the Lemma. 



^ E El- 
^ E E I 



(189) 



□ 



NOTATIONS 



{x}= {x, eZ^lj e{l,...,n}} 

L 

Al = [-L,L]r\Z'^ 
G(")(f}) = {x= (a;i,...,x„)|x, e r!,Vj} 
C'"'(17, m) = {x e G("'(17) I = w for some j} 
C^"\n, 5) = {x e G(")(17) I e for some j} 
d'-^^M = {w e Ml min |g - wl = 1} 

min |g — wj < m — 1} 



aMil/ = {r = (n, ...,r„) € C("'(A/) 3s G Z("'"\G("' (A/) such that max{|ls, -rj\\i<m\ i,j € {1, 

U {r = (ri,...,r„) e Z("'')\G(")(^'^) 3s e C^"^M) such that max{\\si - rj\\i < m\i, j e {1, ...,n}}^ 



diam(x) = max — Xh\ 

j.ke{l,...,n}' 

Ci^^Hn) = {x e G(")(r!) I diam(x) < r} 

Z(x) = max min \x^ — Xk\ 

J,K jeJ,keK 

JUK={l,...,n} 



dist^^(x, y) = max < max dist(xi, {y}), max dist({x},yi) 

\l<i<k l<i<k 

dist^''^^(x,y) = max{dist-H(x,7,y,/),dist?i(xK,yx)} 
maxsuppiJo = [J [suppi?o(k, •) ~ 



EH|Go(x,y)n = ^^E [|G[r^(x,y;i?)| 



dE 
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EXPLANATION OF THE SIMULATION PROCEDURE 



111 this section, we describe the numerics behind the transition from the superfluid phase to the Bose glass phase 
for non-vanishing defect density. The procedure is a standard finite size scaling study in combination with quantum 
Monte Carlo simulations. 

We consider the system with Hamiltonian (in standard lattice notation) 

H = -tY^blb,-J2if^~e,)n, (190) 

for hard-core bosons with tunneling amplitude t and chemical potential /i. The are iid distributed according to a 
uniform distribution [—A, A]. The unit is the hopping t — 1. 

For densities n = and n = 1 the system is a band insulator (empty or full). In the absence of disorder, the system 
is superfluid for any density at zero temperature. Our goal is to determine the critical disorder Ac as a function of 
the defect density for which the system goes over from the superfluid into an insulating Bose glass phase. The Bose 
glass phase has the counterintuitive properties that it is a compressible, gapless though gapless phase Q- Interest 
is especially towards low densities. The system has a "particle-hole" symmetry, so we only need to look at densities 
n < 0.5. 

We will use path integral Monte Carlo simulations that use the worm algorithm [s^ . here in the implementation 



of Ref. l35|. The transition from the superfluid phase to the Bose glass phase is described by the vanishing of the 

(191) 



superfluid density ps which is related to the winding number through the formula [45 



dp 

This reduces to = {W'^)/{2P) in 2d. Note that we know from a general analytical argument that the Bose glass 
always intervenes between the superfluid phase and the insulating vacuum or band insulator [sij . The phase diagram 
of the above Hamiltonian has recently been mapped out for density n = 1 in 3d for soft-core bosons [46[. In 2d, the 
phase diagram has not been published yet. 

The quantum phase transition can be determined numerically from a proper finite size scaling analysis by letting the 
system size L — > oo and temperature T — > 0. From Ref. Q we know that the dynamical critical exponent z = d ~ 2, 
based on the finite value of the compressibility. We can then study the quantum phase transition by taking different 
system sizes L and scaling j3 ^ L'^. We also know the following scaling behavior of the superfluid density in the 
vicinity of the critical point: 

Ps = C\fs{ilL) = L-'U6L'^n, (192) 

because ps ^ ■ Here, \6\ denotes the dimcnsionlcss detuning from the critical point, ^ denotes the correlation 
length, fs is a universal scaling function, and v is the critical exponent. The second equality follows from the fact 
that the correlation length is cut off by the system size if it exceeds the system size. If we scale (3 ^ LF' then the 
curves PsL/^ for different system sizes should intersect in a single point, provided we are in the scaling regime and 
that the irrelevant terms are sufficiently weak. It is computationally advantageous to scale /3 ~ L. Then, the curves 
{W'^) ^ psP ~ psL will not intersect in a single point. Next, we determine the intersection points between curves for 
consecutive system sizes, and extrapolate those intersection points to infinity. Both methods should of course lead to 
the same critical point. 

We have studied the quantum phase transition at fixed density as follows. For every disorder realization we 
determined the chemical potential such that the average number of particles corresponds to our target density. We 
then computed the superfluid density as a function of the bound A for different system sizes L = 4, 8, 16, 24, 32. We 
typically averaged over 100 disorder realizations for the largest system sizes. In the first set of simulations, we scaled 
the inverse temperature as /3 = L^/16, in the second set of simulations we scaled /3 = L, such that the simulations are 
identical for L = 16. The second set of simulations was needed because we did not find a nice single intersection point 
in the first set of simulations in general, hinting that our temperatures are too high or that there may be strong finite 
size effects, which were difficult to overcome with /3 ^ L'^. In these simulations it is computationally better to use the 
scaling /? = L, as well as to have a check on the data and the extrapolations. We then look for the intersection points 
between consecutive curves. Representative data are shown in Fig. [3] for n = 0.125 for quadratic scaling (z = 2) and 
in Fig. IHfor linear scaling (z = 1). The intersection points are then shown in Fig. [6] from which the quantum phase 
transition point can be determined quite accurately as Ac/t = 6.3(2). By repeating this procedure for the densities 
n = 1/100, n = 1/16, n=l/4,n = 3/8 and n = 0.5 we arrive at the phase diagram shown in Fig. [T] 
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A/t 



FIG. 3: (Color online). Curve for the square of winding number as a function of the disorder bound for n = 0.125. Inverse 
temperature is scaled as /3 = 




23456789 
A/t 

FIG. 4: (Color online). Curve for the square of the winding number as a function of the disorder bound for n = 0.125. Inverse 
temperature is scaled as /? = L. 




5.5 6 6.5 7 

A/t 

FIG. 5: (Color online). Same as in Fig. |4l but zooming in on the relevant intersection points. 
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5 I , , , , , 1 

0.05 0.1 0.15 0.2 0.25 0.3 
1/L 

FIG. 6: (Color online). Intersection points for n = 0.125. From these curves, the thermodynamic critical point is found to be 
Ac/t = 6.3(2). The data for z = 1 seem compatible with linear extrapolation, except for the smallest system size which is not 
in the scaling regime. The data for z = 2 are compatible with this extrapolation. 




0.2 0.4 0.6 0.8 1 



n 

FIG. 7: (Color online). Shown is the transition from the superfluid (SF) to the Bose glass (BG) as a function of the density n 
for a system of hard-core bosons on a square lattice and no further interactions between them at T = in the thermodynamic 
limit. 
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